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We calculate the B-parameters for operators arising in theories of new physics beyond the standard 
model (BSM) using HYP-smeared improved staggered fermions on the MILC asqtad lattices with 
Nf = 2 + 1 flavors. We use three different lattice spacings (a w 0.045, 0.06 and 0.09 fm) at a fixed 
ratio of light to strange quarks, me jm s = 1 /5, to obtain the continuum results. Operator matching 
is done using perturbative matching at one-loop order, and results are run to 2 or 3 GeV using 
two-loop running in the MS scheme. For the chiral and continuum extrapolations, we use SU(2) 
staggered chiral perturbation theory. We present preliminary results with only statistical errors. 
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1. Introduction 

In the standard model, mixing in the neutral kaon system arises due to the weak interaction. 
Integrating out the heavy particles, the mixing is described by the matrix element of a AS = 2 
four-fermion operator (Q^ ont below), and is parametrized by Bk- Bk is now determined with high 
precision from lattice QCD [jlj |2^]3[ Q], and plays an important role in constraining the parameters 
of the CKM matrix. In BSM theories, additional operators contribute to kaon mixing. If the matrix 
elements of these operators were known, one could constrain the parameters of these theories in 
a way that is complementary to direct searches. Here we present a calculation of the new matrix 
elements using HYP-smeared staggered valence fermions on the MILC asqtad lattices. 

We adopt the operator basis used in perturbative calculations of anomalous dimensions [||] 

Q C i oaX = [^%(l-75)d a ][s b 7,(\-7 5 )d b ], (1.1) 

^ont = _ y5 )^][^(l _ 75 )^], (1.2) 

^ont = ^^ y(1 _ y^a^^ _ y5)d b^ (1 3) 

2 Cont = _ 75 )^ (1 + (L4) 

Qs° nt = [^7m(1 -Y5)d a ][s b Y^ + Y5)d b ], (1.5) 

where a^ v = [7/i)7v]/2 and a, b are color indices. Q\ leads to Bk, while Q2-5 are the BSM 
operators. The corresponding BSM B -parameters are defined as 

Bi = NIY Vl,LlL-!^S N ^ N *> N5 ^ = ^/ 3 ' 4 ' ~ 2 ' • (L6) 

Ni{K \sY 5 d\0) {0\sY5d\K Q ) 

We stress that this basis of operators differs slightly from the "SUSY basis" used in other lattice 
calculations [§, 0]. We prefer the basis of Ref. [||] since we use perturbative matching and running. 

2. Methodology and Results 

We use the MILC lattices listed in Table |, setting the scale using r x = 0.31 17(6)(+3j) fm 



[gfl. For the valence quarks, we use HYP-smeared staggered quarks [10|, with parameters chosen 
to remove Gicr) taste-symmetry breaking at tree level. Our valence d and s quarks have masses 
denoted m x and m y , respectively, for which we use 10 different values, 

am xy =am s xn/ 10 with n = 1,2,3, 10, (2.1) 



Table 1: MILC lattices used here [|8j]. a is the nominal value of the lattice spacing, "ens" and "meas" are the 
number of gauge configurations measurements per configuration, respectively. ID is an identification tag. 



a (fm) 


ami/am s 


size 


ens x meas 


ID 


0.12 


0.01/0.05 


20 3 x 64 


671 x 9 


C3 


0.09 


0.0062/0.031 


28 3 x 96 


995 x9 


Fl 


0.06 


0.0036/0.018 


48 3 x 144 


749x9 


SI 


0.045 


0.0028/0.014 


64 3 x 192 


747 x 1 


Ul 
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Figure 1: B2(M = 1 l a ) as a function of T = t — 1\. (Red) crosses are from the coarse ensemble C3, with 
(am x ,am y ) = (0.005,0.05); (blue) diamonds are from the fine ensemble Fl, with (am Xl am y ) = (0.003,0.03); 
(purple) octagons are from the superfine ensemble SI, with (am x ,am y ) — (0.0018,0.018); and (brown) 
squares are from the ultrafine ensemble Ul, with (am x ,am y ) = (0.0014,0.014). 

with m s is the nominal sea strange quark mass given in Table |I[ 

The methodology of the calculation for the BSM B -parameters is very similar to that used 
for Bk JTT]] . Many details of the lattice operators and the perturbative matching are given in 
Refs. [Jl3p, although some additional subtleties related to the use of the new operator basis have 



led to small changes Jl 4Q . These, together with the renormalization group running, will be ex- 
plained in Ref. [15]. The kaon and anti-kaon are produced using U(l)-noise wall-sources placed at 
timeslices t\ and ti>t\, while the four-quark operators (and bilinears needed for the B-parameters) 
are placed at an intermediate time t. The resulting B-parameters should be independent of t when t 
is far enough from the wall-sources, so that contamination from excited states is small. Hence we 
fit the data to a constant in the plateau region. The fitting range is determined using the two-point 
correlator from the wall-sources to the taste-^ axial current. In Fig. [j], we show results for B 2 
as a function of T = t — 1\ with our most physical kaon. When fitting, we ignore the correlations 
between timeslices (diagonal approximation for the covariance matrix) to avoid an instability of 
the fit due to small eigenvalues of the covariance matrix. The fitting errors are estimated using the 
jackknife method. 

To increase statistics, we perform multiple measurements with randomly chosen t\ on each 
gauge configuration (see Table [I]). We find considerable autocorrelation for the BSM B-parameters 
on the fine, superfine, and ultrafine ensembles. Hence we bin the data on these ensembles, using a 
bin size of 5. 

After calculating the BSM B-parameters for 55 valence quark mass combinations, we perform 
the chiral extrapolation to the physical down and strange quark masses. We first extrapolate m x to 
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Figure 2: (a) = 1 AO from the NNNLO Bayesian X-fit vs. X P , on Fl, for am y = 0.03. The red diamond 
represents the physical point, (b) B2(n = 1 /a) from the Y-fit vs. Yp, on Fl. The red diamond corresponds to 
the physical strange quark mass. 



m^ hys at fixed m y ("X-fit"), and then linearly extrapolate m y to m? hys ("Y-fit"). In the X-fit, we fit 
to the form from SU(2) staggered chiral perturbation theory (SChPT), which requires m x <C m y . 
Hence we take lightest four quark masses for m x (e.g. m x = {0.005, 0.01, 0.015, 0.02} on the 
coarse ensemble) and the heaviest three quark masses for the m y (e.g. m y = {0.04, 0.045, 0.05} 
on the coarse ensemble). 

For the X-fit we use the next-to-leading order (NLO) SChPT result for the Bj from Ref. 
extended to higher order: 

fl^X-fit, NNNLO) = ciFoO") + c 2 X + c 3 X 2 + c 4 X 2 ( ln(X)) 2 + c 5 X 2 ln(X) + c 6 X 3 . (2.2) 



Here X = Xp/At, with Xp the squared mass (in physical units) of the taste-^5 pion composed of 
two light quarks, Xp = M 2 X . P . For the chiral renormalization scale we take A x = lGeV. Fo(j) 
contains the leading order and NLO chiral logarithms, and is completely known in terms of f n 



and measured lattice pion masses [13]. The c 2 term is the NLO analytic term. We include three 



"generic" NNLO terms: the C3 term whis is representative of NNLO analytic terms, and the C4 and 
C5 terms, which are representative of NNLO chiral logarithms in continuum ChPT. We also include 
one NNNLO term, with coefficient c&. 



We fit using the Bayesian method Jlq ] with parameters C4_6 constrained to be of order unity, 
which is the expectation from chiral power-counting. Specifically, we constrain them to be £4-6 = 
Oil. The full correlation matrix is included in the X-fit. 

Having determined the parameters ci_6, we can simultaneously extrapolate the results to the 
physical point m x = m^ hys and remove lattice artifacts due to taste-breaking in pion masses in the 
chiral logarithms Fo(j), as explained in Ref. [11]. We also set mi 7^ m? hys in the logarithms. In 



Fig. 2(a), we show the X-fit for B 2 on ensemble Fl. 
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After the X-fit, we perform the Y-fit, in which we extrapolate m y to the physical strange quark 
mass mj? ys . We expect that the Bj are smooth, analytic functions of Yp. It turns out that a linear 
form describes the data well: 



Bj(Y-&t) = b l +b 2 Y P , 



(2.3) 



where Yp = M 2 vP is the squared mass of the valence pion with composition yy and taste ^5. Fig. 2(b) 
shows results of the Y-fit for B2 on Fl. At this stage, we use uncorrected fitting for the Y-fit. 

After the chiral extrapolations, we know the BSM B-parameters evaluated at a fixed lattice 
spacing and matched to the MS scheme at a scale pL = \/a. In order to extrapolate to the continuum 
limit {a = 0), we need to first run the results to a common scale ju. In the RG running, operator 
mixing arises in pairs: {Q% oat , Q^ ont ) and (Q% ont , Q^ ont ). 

The anomalous dimension matrix for the BSM AS = 2 operators in the basis of Eqs. (L2)-(L5) 
is calculated up to two-loop order in Ref. This is in the MS scheme with naive-dimensional 
regularization of 75 and with the choice of evanescent operators made by Ref. f^]. 1 Hence we 
calculate the RG evolution matrix for the BSM B-parameters at that order. In the case of RG 
running for B4 5 , there is a removable singularity in the standard two-loop approximate analytic 
solution. To resolve this, we use the analytic continuation method introduced in Ref. Ji7|]. We have 
checked the results by numerical evolution of the RG equations. The resulting BSM B-parameters 
evaluated at jU = 2 GeV and 3 GeV are given in the Tables [2| and ||. We note that statistical errors in 
the BSM B-parameters are smaller than those in Bk- 

The final step is the continuum extrapolation of the results. We know that the leading a and GC S 



dependence to be [14] 



Bj = di+d 2 (aA) 2 + d 3 (aA) 2 a s + d 4 a 2 + d 5 (aA) 4 + • • 



(2.4) 



where a s = a, , (1/a). We do a Bayesian fit to this form, taking the QCD scale determining the 
magnitude of discretization errors to be A = 300MeV, and constraining d 2 -5 = ± 2. As for Bk, 
we find that fits to all four lattice spacings are very poor, with ^ ug /dof = 6.6 ~ 30 for #25- Thus 
we drop the results from the coarse lattice and fit to the finest three spacings (Fl, SI and Ul). In 



Table 2: Preliminary results for BSM B-parameters and Bk at jJ. 
using linear extrapolation. Only statistical errors are shown. 



2 GeV. Continuum values are obtained 



B 7 \Lat 


C3 


Fl 


SI 


Ul 


Continuum 


B K 


0.5672(52) 


0.5295(43) 


0.5362(38) 


0.5318(70) 


0.5383(66) 


B 2 


0.5404(09) 


0.5646(14) 


0.5967(19) 


0.6058(31) 


0.6245(30) 


B 3 


0.3689(06) 


0.4148(10) 


0.4594(14) 


0.4805(24) 


0.5032(22) 


B 4 


1.0965(23) 


1.1260(28) 


1.0911(37) 


1.0942(57) 


1.0698(56) 


B 5 


0.9278(20) 


0.9381(25) 


0.8875(31) 


0.8720(49) 


0.8432(48) 



'A different choice of evanescent operators was made in the one-loop matching calculation of Ref. Jl^]. We have 
now extended this calculation to the scheme of Ref. [ph. Results will be reported in Ref. [13]. 
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Table 3: Preliminary results for BSM B-parameters and Bk at fi = 3GeV. Notation as in Table ^[ 



B 7 \Lat 


C3 


Fl 


SI 


Ul 


Continuum 


B K 


0.5478(50) 


0.5114(42) 


0.5179(37) 


0.5137(67) 


0.5199(64) 


B 2 


0.4779(08) 


0.4993(12) 


0.5277(17) 


0.5358(28) 


0.5524(26) 


B 3 


0.3152(05) 


0.3496(08) 


0.3840(12) 


0.3997(20) 


0.4174(19) 


B 4 


1.0462(22) 


1.0750(26) 


1.0421(36) 


1.0452(55) 


1.0222(54) 


B 5 


0.9132(19) 


0.9272(24) 


0.8824(31) 


0.8714(48) 


0.8450(47) 




Figure 3: Continuum extrapolation of B2 at 2GeV. (Red) diamond is the result from the linear fitting func- 
tion; and (blue) circle is the result from the Bayesian constrained fitting with the fitting function given in 
Eq. @. 

this case, both the linear fitting (keeping only d\ and d-i) and the constrained fitting (with J1-5) 
work well. In Fig. [§, we show an example of the continuum extrapolation for #2- In Tables || and 
|H we quote the results from the linear extrapolation. Clearly the systematic errors associated with 
the choice of continuum extrapolation are significantly larger than the statistical errors. 

3. Outlook 

The next stage in our calculation is to quantify all sources of systematic error and so draw up a 
complete error budget. This requires results at other values of the light sea-quark masses to estimate 
residual ma dependence, and at other volumes to estimate finite volume effects. The latter can also 
be estimated using SChPT, and are expected to be small. We also plan to investigate whether the use 
of ratios which cancel chiral logarithms reduces errors in the analysis, and to compare the results 
to those from an analysis using SU(3) SChPT. We expect that, as for Bk, our dominant errors will 
come from to the use of one-loop matching and the continuum extrapolation. 
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Although our results are preliminary, it is interesting to compare them to those found in 
Refs. [@, 0]. Changing to the SUSY basis, B 2 , B 4 and B 5 are unchanged, while fif USY = (5B 2 - 
3Bt,)/2. Thus our preliminary results in the tables translate into B^ VSY = 0.81 and 0.75 at pL = 2 
and 3 GeV, respectively. There are some disagreements between our results and those of Refs. ^ 0] 
at the 25% level. Determining whether these are significant will require our full error budget. 
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